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ABSTRACT 

Many astrophysical bodies harbor magnetic fields that are thought to be sustained 
by a dynamo process. However, it has been argued that the production of large-scale 
magnetic fields by mean-field dynamo action is strongly suppressed at large magnetic 
Reynolds numbers owing to the conservation of magnetic helicity. This phenomenon is 
known as catastrophic quenching. Advection of magnetic fields by stellar and galactic 
winds toward the outer boundaries and away from the dynamo is expected to alleviate 
such quenching. Here we explore the relative roles played by advective and turbulent 
diffusive fluxes of magnetic helicity in the dynamo. In particular, we study how the 
dynamo is affected by advection. We do this by performing direct numerical simula- 
tions of a turbulent dynamo of a 2 type driven by forced turbulence in a Cartesian 
domain in the presence of a flow away from the equator where helicity changes sign. 
Our results indicate that in the presence of advection, the dynamo, otherwise station- 
ary, becomes oscillatory. We confirm an earlier result for turbulent -diffusive magnetic 
helicity fluxes that for small magnetic Reynolds numbers (Rm < 100. ..200, based on 
the wavenumber of the energy-carrying eddies) the magnetic helicity flux scales less 
strongly with magnetic Reynolds number (Rm -1 / 2 ) than the term describing magnetic 
helicity destruction by resistivity (Rm - ) . Our new results now suggest that for larger 
Rm the former becomes approximately independent of Rm, while the latter falls off 
more slowly. We show for the first time that both for weak and stronger winds, the 
magnetic helicity flux term becomes comparable to the resistive term for Rm > 1000, 
which is necessary for alleviating catastrophic quenching. 
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1 INTRODUCTION 

A theoretical framework for explaining the large-scale mag- 
netic fields observed in many astrophysical bodies is mean- 
field dynamo theory. Its basic idea is that the inductive ef- 
fects of turbulent motions are able to amplify a weak mag- 
netic field and maintain it on timesca les longer than the 
magnetic diffusion time (|Moffatdll978h . Gradients in the 
large-scale velocity field, like shear motions, can also con- 
tribute significantly to the amplification of the magnetic 
field. In mean-field dynamo theory the contribution of the 
turbulent scales is parameterized through the electromotive 
force which depen ds on the large-scale ma gnetic field as well 
as its derivatives (|Krause fc Radlerlll~980iy The coefficients 
in front of the magnetic field and its derivatives are called 
turbulent transport coefficients. They can describe either 



* E-mail: fabio@nordita.org (FDS) 
© 2002 RAS 



turbulent-diffusive (with turbulent diffusion cc %) or non- 
diffusive (e.g., the a effect or turbulent pumping) effects. 

Under some approximations (e.g., in the low conductiv- 
ity limit for small magnetic Reynolds number, Rm 1, or in 
the high conductivity limit for short correlation times, i.e., 
small Strouhal number, St ^ 1), theories like the first order 
smoothing approximation are able to predict the functional 
form of the expressions and the correct values of the coeffi- 
cients. Within their limits of validity, these results present 
a remarkably good agreement with the computation of the 
transport coeffici ents through dir ect numerical simulations 
(DNS); see, e.g., ISur et all (|2008l ). However, not enough is 
known about the functional form of these coefficients at 
large values of Rm (i.e., small values of the microphysical 
magnetic diffusivity) and about the saturation process when 
the magnetic field becomes dynamically important. Under- 
standing the behavior of the dynamo in these regimes has re- 
mained an important problem for several decades. Although 
many recent works have contributed to understanding dy- 
namo saturation at large magnetic Reynolds numbers, more 



2 F. Del Sordo, G. Guerrero and A. Brandenburg 



work is still necessary to have a complete picture of the dy- 
namo excitation and saturation mechanisms. 

Among the turbulent transport coefficients the a ef- 
fect is particularly important, because it allows a closed 
dynamo loop for regenerating both poloidal and toroidal 
magnetic fields. It has been suspected, however, that in 
closed or triply periodic domains the a effect can be 
strongly suppressed at higher mag netic Reynolds numbers 
and might scale like a oc Rm " 1 (|Vainshtein fc Cattaned 
ll992l ; ICattaneo fc H ughes 199 6]). An expla nation for this was 



proposed by Gruzinov fc Diamond! (I1994T ). who used the a 



effect derived by Pouquet et al.l ( 19761 ). which has, in ad 



dition to the kinetic helicity density, a contribution propor- 
tional to the current helicity density. It is this quantity which 
builds up as the dynamo saturates. 

This is a consequence of magnetic helicity conservation 
and can be explained as follows: the large-scale magnetic 
field generated by the a effect is helical, but in order to sat- 
isfy the conservation of total magnetic helicity, a small-scale 
field with equally strong magnetic helicity of opposite sign 
must be generated in the system. The small-scale magnetic 
helicity is responsible for the creation of a magnetic a effect 
(oh) which contribu tes with opposite sign to the kinetic 
a. This basic idea led iKleeorin fc Ruzmaikm] (|l982h to pro- 
pose the dynamical quenching model at a time well before 
simulations saw any indications of catastrophic quenching. 
Even nowadays the issue is quite unclear when it comes to 
making predictions about the high-Rm regime. The final 
amplitude that the magnetic a effect acquired depends on 
the geometry of the system and on the value of the mag- 
netic Reynolds number. For highly turbulent astrophysical 
objects with high Rm like the Sun or the Galaxy, «m could 
attain higher amplitudes, decreasing then the dynamo effi- 
ciency. However, the dynamics of «m also depends on the 
ability of the system to get rid of the small-scale magnetic 
helicity responsible for its creation. In a closed or triply pe- 
riodic homogeneous domain, magnetic helicity annihilation 
depends just on the microscopic magnetic diffusivity. This 
is a very slow process given the scales and diffusivity val- 
ues under consideration. However, an obvious solution to 
this catastrophic (Rm- dependent) quenching is to allow the 
system to get rid of helical small-scale magnetic fields. 

In real astrophysical systems, this processes is generally 
expected to happen in a number of different ways. Among 
the various mechanisms for removing magnetic helicity from 
the system we focus here on the role played by the turbulent- 
diffusive magnetic helicity flux and by the presence of advec- 
tive flows or winds. The role of these magnetic helicity fluxes 
has been tested in the context of mean-field dynamo mod- 
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models have demonstrated the importance of magnetic he- 
licity fluxes in solving the catastrophic quenching problem. 

Verifying the validity of these results in DNS is more 
complicated since obtaining higher Rm in the numerical 
models requires high resolution and large computational 
resources. Various attempts have, however, succeeded in 
demonstrating the role of magnetic helicity c onservation in 
the sa turation of the dynamo. For instance, iBrandenburd 
i|200ll) studied the saturation in triply periodic helically 



forced dynamos of a 2 type. The role of open magnetic 
boun dary condition s for convect i ve dy namos has been stud- 
ied in lKapyla et alj (120081. l200gj. l2010h . Furthermore, using 
forced turbulence. iMitra et all <|2010ah (hereafter MCCTB) 
have verified the existence of turbulent-diffusive magnetic 
helicity flux es in a 2 dynamo models i n the presence of an 
equator and lHubbard fc Brandenburd (|2010h (hereafter HB) 
did the same for a dynamo region embedded inside a highly 
conducting halo which provided a more realistic boundary 
condition. In both cases it was found that a fit to a Fickian 
diffusion law can account for this flux and that the diffusiv- 
ity value is comparable to or below the value of the turbulent 
magnetic diffusivity %. The resulting Fickian diffusion coef- 
ficient was found to be approximately independent of Rm. 
By considering a statistically steady state, and noting that 
the local value of the magnetic helicity density was also sta- 
tistically steady, their result became then also independent 
of the gauge chosen to define the magnetic vector potential. 
In addition, shear flows have been argued to be effec- 



tiyein alleviating catastrophic quenching (jVishniac fc Choi 
l200lf ) and allowing signifi cant saturation levels of the dy- 
namo l|Kapvla et alJuOOSl ). although it appears now plau- 
sible that their result could also be explained through a 
change in the excitation conditions of the dynamo. Indeed, 
recent DNS have f ailed to demonstrate the presen ce of the 
Vishniac-Cho flux (|Hubbard fc Brandenburdl201ll ). Yet an- 
other possibility is the advective magnetic helicity flux. In 
the context of the galactic dynamo, alleviation of catas- 
trophic quenching th anks to a wind has b een studied in 
mean- fi eld models bylShukurov et alj (|2006h and | ur et al.l 
(|2007h . IMitra et all (|201lh studied the role of a wind in so- 
lar mean-field dynamo models. The models studied in the 
present paper allow us to compare with their results and to 
determine the importance of magnetic helicity fluxes in the 
dynamical evolution of the magnetic a-effect. To our under- 
standing the study of advective fluxes in DNS of a dynamo 
is an outstanding problem. With this paper we intend to 
close this gap. 

We perform DNS leading to a 2 -type dynamo action in a 
domain with kinetic helicity of opposite signs on both sides 
of the equator. We use a relaxation term to include a large- 
scale flow that advects the large-scale magnetic field. Fur- 
thermore, we consider periodic boundary conditions in the 
horizontal directions, zero-gradient conditions for the veloc- 
ity and vertical field conditions for the magnetic field. In this 
way we allow for the removal of magnetic helicity through 
advection. For the sake of simplicity and to study these ef- 
fects separately in a clear way, we do not include large-scale 
shear. Nevertheless the results presented here should also be 
applicable in the context of the galactic dynamo and, in prin- 
ciple, also to the solar dynamo, where large-scale winds have 
been shown in mean-field models to play a role in carr ying 
magnetic helicity outside its bounds (|Mitra et alj|201ll ). 

This paper is organized as follows. In Sect. [2] we de- 
scribe the physical model considered here and present the 
equations governing its evolution. In Sect. [3] we present the 
results of the simulations. First we describe the properties 
of the solutions without the wind. Next, we explore the ef- 
fects that the wind has on the characteristics of the dynamo 
solution. Finally, we determine the magnetic helicity fluxes 
present in the model and verify their balance with the pro- 
duction terms to prevent the quenching of what corresponds 
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to the a effect in the related mean-field description. We con- 
clude and summarize the results in Sect. [4] 



2 THE MODEL 

2.1 Governing equations 

We use the Pencil Code0 to solve the following set of com- 



pressible hydromagnetic equations in an isothermal layer 
OA 
~dt 



U x B — fj,orjJ, 



(1) 



^ = - C 2 S V hip + - J x B + i V • 2upS + f w + f, (2) 
Dt p p 



Dlnp 
Dt 



= -V -U + q p , 



(3) 



where D/Dt = d/dt + U ■ V is the advective derivative, A 
is the magnetic vector potential, B — V x A is the mag- 
netic field, J — V x B/p,o is the current density, /io is the 
magnetic permeability, r\ and v are magnetic diffusivity and 
kinematic viscosity, respectively, c a = const is the sound 
speed, U is the velocity, p is the density, S is the rate of 
strain tensor given by 



s« = §(Ui,i + tfj,i)-!*«v-tf ) 



(4) 



where the commas denote derivatives, / w provides a forcing 
for the wind (defined below in Section I2.2|l , q p is a source 
term in Eq. ([3| needed to replenish the resulting mass loss, 
/ is a time-dependent random <5-correlated forcing function 
of the form 



f = f(x,t;a(z)), 



(5) 



where a is related to its local helicity density, 

(f-Vxf)/{k l f) = 2a/{l + a 2 ), (6) 

and is chosen to vary like cr(z) = sin(2n z/L z ) with a 
sign change across the equator at z = 0. This forcing 
drives turbulence in a band of wavenumbers around fcf. The 
modulation a(z ) of th is forcing is similar to that used by 
IWarnecke et al.l (|201ll ) to simulate a sign change of helicity 
in forced turbulence in a spherical wedge. 

We consider a computational domain of size L x x L y x 
L z , with quadratic horizontal extent, L x = L y , using peri- 
odic boundary conditions and a vertical extent that is twice 
as big, L z = 2L X , with \z\ L z /2 (i.e., -L z /2 < z < L z /2) 
and an equator at z — 0. Our boundary conditions are 



Ux , Z Uy j z Uz 



t/w = A X ,z = Ay,z = Az = (7) 

on the top and bottom boundaries at z = ±L./2 = ±z top - 
(7 W is the wind profile, defined below in Section 12.21 The 
lowest horizontal wavenumber in the domain is k\ — 2tt/L x . 
In the following, we use ki as our inverse length unit, so 
\kxz\ ^ 27r. To eliminate boundary effects, we restrict most 
of the analysis to a diagnostic layer, \z\ ^ L* with fciL, = 3. 
For all our runs we choose fcf/fci = 4, which is a compromise 
between it being large enough to allow a large-scale mag- 
netic field to be generated and yet small enough to achieve 
sufficiently large values of Rm. 
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We set c s to unity in the code, so our dimensionless time 
is in units of the sound travel time, (c s fci) _1 . However, the 
relevant physics is not governed by compressibility effects, 
so it is more natural to quote time in turnover times, i.e., 
we quote instead the value of tu TUls k{. In most of the cases 
reported below, the turbulent Mach number, Ma = u rms /c s , 
is around 0.1. Likewise, in the code v and r\ are given in 
units of Ca/ki, but it is physically more meaningful to quote 
corresponding Reynolds numbers instead. Our resolution is 
increased with increasing values of Rm, so the largest reso- 
lution used in this paper is 1024 x 1024 x 2048 meshpoints. 
We return to this issue at the end of the paper. 

2.2 Generating the wind 

In our model, the advective term from the wind is given by 
the forcing function in Eq. |(5J, 

U = -—\U-U*(zj\, (8) 

Tw 

where U is the horizontally averaged velocity field, and 



U v (z) = U - 



•2-top 



(9) 



is the wind profile that increases linearly toward the z 
boundaries. The wind profile can be modified by the turbu- 
lence and the magnetic field, but the original outflow profile 
tends to be restored on a timescale r w . The presence of a 
wind leads to mass loss across the vertical boundaries with 
a mass loss rate that depends on Uq. 

Stellar winds are the main agents of mass loss in stars. 
In a galactic environment it is possible to observe galactic 
winds as well as galactic fountains. These mechanisms can 
be driven by the explosions of supernovae in the galactic 
disc. In this case a direct estimate of the mass loss rate 
is more complicated, given that it is expected to be very 
small. However, to have stationary conditions, we keep the 
mass in the domain constant using the source term q p in 
Eq. ©. This source term tends to be restored the density 
at each spatial point in the domain to its initial value po on 
a timescale r s = r w . Thus, analogously to Eq. (JS|, we write 
q p = -r^Qxip-lnpo). 

We study the dependence of our model on the dimen- 
sionless wind speed and the magnetic Reynolds number of 
the turbulence. These are defined as 

c ^ ' ^ w T)™ u rms /ir\\ 

bw = t— , Rm = — — . (10) 

In all cases, we use a magnetic Prandtl number of unity, i.e., 
v/r) = 1. 

2.3 Magnetic helicity fluxes 

In our model we expect two different kinds of magnetic he- 
licity fluxes: those caused by the wind, i.e. advective mag- 
netic helicity fluxes, and those due to turbulence in the pres- 
ence of a mean gradient of the magnetic helicity density, i.e. 
turbulent-diffusive magnetic helicity fluxes. To assess their 
importance in the magnetic helicity budget, we now consider 
the magnetic helicity equation in the Weyl gauge which is 
used in Eq. JT]), i.e., 

^A~B = -2 W0 JVB- V-^, (11) 
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where overbars denote averages over x and y and J- = 
£ x A is the total magnetic helicity flux, with E = 
77/10 J U x B being the electric field in the lab frame. 
This equation is evident ly gauge-dependent; see for instance 
ICandelaresi et all (|201ll ). In particular, since A ■ B is not a 
physical qu antity, it could drift - eve n in the steady state; 
see Fig. 2 o f lBrandenburg et al.l (|2002h for an example. How- 
ever, if A ■ B is constant in a particular gauge, then we have 

V ■ 7 = -2rm ~J~B, (12) 

where now V- J- must be gauge- independent, because J and 
B are gauge-invariant. This argument was invoked by MC- 
CTB and HB to determine turbulent-diffusive contributions 
to the magnetic helicity flux. 

In the present work, we are interested in two contribu- 
tions to h = A ■ B, one from the mean fields, h m = A ■ B, 
and one from the fluctuating fields, hi = a ■ b. Their sum 
gives the total mean magnetic helicity density, i.e., h = 
hm + hf. Note, however, that only hf is the component di- 
rectly relevant for the study of catastrophic quenching, be- 
cause it is approximately proportional to the current helicity 
density, j ■ b, which in turn determines the magnetic contri- 
bution to the a effect. [The approximate proportionality of 
magnetic and current helicities is non-trivial and will need to 
be re-assessed below; see also Fig. 3 of MCCTB and Table 2 
of HB for earlier examples.] 

The evolution equation for hi is 

^ = -2S-B-2 Vl i j~b-V-7 f , (13) 

where, as mentioned above, we allow two contributions to 
the flux of magnetic helicity from the fluctuating field tFf. 
an advective flux due to the wind, J-i" = hfU w , and a 
turbulent -diffusive flux due to turbulence, modelled here 
by a Fickian diffusion term down the gradient of hi, i.e., 
jFf dlff = — KhV/if. Here, £ = u x b is the electromotive 
force of the fluctuating field. 

In the steady state, and if h{ is then also constant (which 
is not guaranteed to be the case because hi is a priori gauge- 
dependent), we have 

V • J^t = —2£ • B — 2j]fM) J~T>- (14) 

Again, although V • J-i is in principle gauge-dependent, it 
can now be determined by measuring £B and j ■ b, that are 
manifestly gauge-independent quantities. This means that 
V ■ Ti must be gauge-independent as well. We assume that 
J-i has a component only in the vertical direction. We can 
therefore obtain its z dependence through integration via 

J fe = / V-Ttdz'. (15) 
Jo 

The assumption of only a z component of J-i would 
break down in the presence of shear, where cross- 
stream fluxes with fini t e div ergence are possible; see 
iHubbard fc Brandenburg! i|201lT ). 

For the discussion of our results presented below, let us 
contrast our present simulations with those of MCCTB. In 
their case, the outer boundary condition at z — ±L z /2 was 
a perfect conductor (P.C.) one and the most easily excited 
mode was antisymmetric about the midplane with dynamo 
waves propagating toward the equator. This antisymmetry 
results in permitting a flux of magnetic helicity through the 



Table 1. Comparison of boundary conditions and other proper- 
ties of the simulations of MCCTB and the present work. 





IVlL^o r rS 


Til) 

riD 


present work 


boundary 


P.C. 


halo 


V.F. 




T U = 


T fz + 






V • Pi ± 


V • 3=i = 


V Ft = 


equator/ 


antisymmetry 


symmetry 


symmetry 


midplane 


(like V.F.) 


(like P.C.) 


(like P.C.) 






7iz = 






V J=i = 


V ■ F t ^ 


V J=i = 



equatorial plane and in this sense has the same effect as the 
vertical field (V.F.) boundary condition. This, together with 
the fact that the magnetic helicity density is antisymmet- 
ric about the equator, is the reason why in their case the 
turbulent-diffusive flux can play a measurable role. How- 
ever, because Tt z has vanishing vertical derivative at the 
equator, the V ■ jFf vanishes there. This is different in the 
model of HB, in which the helicity is arranged to be sym- 
metric about the midplane, which is therefore not an equa- 
tor in the usual sense. Here the field is symmetric about 
the midplane, corresponding thus to a P.C. condition, and 
thus V • J-i 7^ 0. The boundary conditions and their prop- 
erties are summarized in Table [T] for MCCTB and HB and 
compared with those used in the present work. 

Unlike MCCTB, in the present work the V.F. condition 
is applied on the outer boundaries, in which case the most 
easily excited mode is symmetric about the equator with dy- 
namo waves travelling away from the midplane. This is sim- 
ilar to a P.C. condition at the midplane, for which the mag- 
netic helicity flux vanishes. However, because h{ is antisym- 
metric about the equator, it must have a turning point there, 
so its second derivative vanishes and V-^f = 0. The present 
model does not have shear, but the nature of the domi- 
nant mode is similar to early simulation s of dynamos driven 
by th e magneto-rotational instability (jBrandenburg et al.l 
Il995h . 



3 RESULTS 

3.1 Model without advective flux 

We begin by describing the results for a dynamo in the ab- 
sence of an advective flux (Sw = 0). The solution for this 
particular setup is a steady magnetic field mainly concen- 
trated around the equator of the domain, where the mag- 
netic helicity changes its sign. In Figure [1] we show the B y 
and B z components of the magnetic field in the saturated 
phase of a model without wind and Rm = 206 (later referred 
to as Model N3). Note that B x = B y — on the top and 
bottom boundaries, owing to the use of vertical-field bound- 
ary conditions. Both of them, as well as B z , do not show 
any significant temporal change once fe rms has reached its 
saturation value. This can be observed in the top panel of 
Figure [5] where the vertical distribution of B y is depicted 
as a function of time. 

The fact that this model is steady in the absence of 
a wind is surpri sing, because according to linear mean- 
field calculations (|Brandenburg et al.ll2009h it should exhibit 
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B/B. 



B v /B n t =2870 B,/B n 



Figure 1. Visualization of B y (left) and B z (right) on the bor- 
ders of the domain for model N3 in the saturated phase of the 
simulation (t is time in turnover times). 



cyclic behavior with dynamo waves moving away from the 
midplane. This discrepancy could be related to nonlinearity 
or to differences resulting from the use of mean-field theory. 
However, for the different boundary conditions used by MC- 
CTB, mean-field and direct numerical simulations exhibit 
rather similar behavior. If it is a consequence of nonlinear- 
ity, it could be related to not allowing magnetic helicity to 
escape the domain. Indeed, the behavior is certainly quite 
different from the cases with advective magnetic helicity flux 
(see below), and it is also different from the otherwise similar 
accretion disc models. 



3.2 Model with advective flux 

Let us now turn to models in which a wind is included 
(Sw 7^ 0). An example of the resulting wind profile as well 
as the vertical distribution of u rms is shown in Figure [3] Even 
with just a weak wind the dynamo becomes oscillatory; see 
Figure [2] Note that the cycle period decreases as the wind 
speed is increased. We observe oscillatory solutions of even 
parity, that is B x and B y are on average symmetric with 
respect to the midplane z = 0, with dynamo waves migrat- 
ing away from z = 0. T his is expected based on m ean-field 
models in similar setups (|Brandenburg et alj|2009l ) provided 
the outer boundary condition is a vacuum or vertical field 
condition, as is the case here. 

In Figure [4] we can see how the actual B y (x,y, z,t), 
as opposed to its horizontal average B y (z,t), evolves dur- 
ing half a period in the saturated phase of the simulations, 
changing gradually from negative to positive polarity. In Ta- 
ble[2]we summarize important output parameters that char- 
acterize the simulations and, in particular, details regarding 
the magnetic helicity balance. Note that all table entries are 
non-dimensionalized by normalizing with relevant quantities 
such as B Bq ; see the table caption for details. Magnetic he- 
licity and the various production terms are antisymmetric 
about the midplane. Within the range \z\ ^ L«, all these 
quantities vary approximately linearly with z. Therefore we 
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Figure 2. Space-time diagrams of B y for different wind intensi- 
ties Sw corresponding to Models N3, W3, M2, and S2 from top 
to bottom. The time axes have been shifted such that for each 
run about 1200 turnover times are being displayed. Note that the 
cycle period decreases with increasing wind speed. 
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Figure 3. Resulting vertical profile of U z together with the 
rms velocity as a function of height. Different lines correspond 
to different times. In this case Uq = 0.015c s , corresponding to 
S w = 0.0055. 



characterize their values by their slope. An appropriate nor- 
malization is therefore kir]toB Bq . 

As can be seen from the bottom panel of Figure [5] the 
difference between the values of total and turbulent-diffusive 
fluxes is roughly constant with z, so that its divergence is 
small. This shows that in this particular setup the turbulent- 
diffusive magnetic helicity flux has actually no contribution 
in balancing the rhs of Eq. (|14|) to zero. This is different 
form the case studied by HB, in which a finite magnetic 
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Figure 4. Visualization of B y (left) at six different times during the evolution of the system, for Model S3. It is visible its variability, 
being this component of the magnetic field prevalently negative in the first snapshot and gradually turning positive in the others. Time 
is given in turnover times and spans over half a cycle. On the right, B z is visualized on the borders of the domain for model S3. It does 
not show any relevant variability during its evolution. 



Table 2. Characteristic output parameters of the simulations. Here, Syy characterizes the wind speed, B rma is the rms value of the 
mean field normalized by B cq , 2& ■ B, 2r)j ■ b, and V • give magnetic hclicity production, dissipation, and flux divergence in units of 
&f?7to5g q , is the turbulent— diffusive magnetic helicity flux in units of r/toB^q, characterized by the diffusion coefficient K{/ijt, fc e ff is 

normalized by fcf , j ■ b is normalized by kfB% q , and V ■ ^F m is the flux divergence of the mean field in units of kfqtoB^. The strongest 
outflows we simulate are those of Models II and 12, for which we reckon that the magnetic field is manly carried out of the domain by 
the outflow. The outflows on which we manly focus in this work are those of Models SI — S6, in which the maximum value of the wind 
speed is Uq « 0.15 ■ it rms . N x indicates the number of mesh points in the x direction. (In all cases we have N v = N x and N z = 2N X .) 



Model 


Sw 


Rm 


2£ ■ B 




2qjb 


V 


Ft 






fceff 


i 


b 




Tl 


0.0000 


9 


0.066 ± 0.019 


-0 


069 ±0.018 


0.004 ±0 


001 


0.007 


-0.3 ± 0.7 


1.22 


-0 


03 


0.06 


T2 


0.0000 


23 


0.032 ± 0.005 


-0 


035 ± 0.003 


0.004 ± 


007 


0.002 


-0.5 ± 0.4 


1.16 


-0 


03 


0.03 


Nl 


0.0000 


37 


0.048 ± 0.007 


-0 


047 ± 0.007 


0.001 ±0 


001 


0.004 


0.1 ± 0.2 


1.48 


-0 


07 


0.05 


N2 


0.0000 


81 


0.023 ± 0.008 


-0 


022 ± 0.005 


0.000 ± 


005 


0.002 


-0.0 ± 0.2 


1.37 


-0 


07 


0.02 


N3 


0.0000 


206 


0.001 ±0.001 


-0 


002 ± 0.000 


0.001 ±0 


001 


-0.003 


0.0 ±0.2 


0.97i 


-0 


01 


0.00 


N 1 


0.0000 


397 


0.000 ± 0.004 


-0 


001 ± 0.001 


-0.000 ±0 


004 


-0.000 


0.1 ± 0.2 


0.84i 


-0 


01 


0.00 


N5 


0.0000 


722 


-0.006 ± 0.002 


-0 


005 ± 0.001 


0.005 ± 


005 


0.017 


0.1 ± 0.3 


2.34 


-0 


16 


-0.02 


N6 


0.0054 


1073 


0.010 ±0.004 


-0 


006 ± 0.000 


-0.018 ±0 


015 


0.019 


0.0 ±0.3 


2.83 


-0 


28 


0.02 


Wl 


0.0020 


24 


0.205 ±0.007 


-0 


196 ± 0.007 


-0.008 ±0 


002 


0.013 


0.4 ±0.1 


1.17 


-0 


19 


0.19 


W2 


0.0019 


51 


0.094 ± 0.022 


-0 


088 ± 0.023 


-0.005 ±0 


001 


0.011 


0.5 ± 0.1 


1.45 


-0 


18 


0.09 


W3 


0.0019 


129 


0.047 ± 0.004 


-0 


043 ± 0.004 


-0.004 ±0 


002 


0.010 


0.4 ±0.2 


1.60 


-0 


23 


0.05 


W4 


0.0018 


265 


0.026 ± 0.002 


-0 


024 ± 0.001 


-0.003 ±0 


003 


0.008 


0.2 ± 0.2 


2.03 


-0 


26 


0.03 


W5 


0.0018 


540 


0.014 ±0.004 


-0 


012 ±0.001 


-0.002 ±0 


012 


0.008 


0.0 ±0.2 


2.62 


-0 


26 


0.01 


M2 


0.0038 


51 


0.090 ± 0.007 


-0 


082 ± 0.010 


-0.006 ± 


001 


0.008 


0.4 ±0.2 


1.48 


-0 


17 


0.09 


SI 


0.0060 


24 


0.167 ±0.004 


-0 


152 ± 0.004 


-0.012 ±0 


002 


0.019 


0.8 ± 0.2 


1.27 


-0 


15 


0.16 


S2 


0.0056 


51 


0.085 ± 0.004 


-0 


074 ± 0.004 


-0.007 ±0 


007 


0.015 


0.5 ± 0.4 


1.52 


-0 


16 


0.08 


S3 


0.0055 


133 


0.034 ± 0.005 


-0 


029 ± 0.004 


-0.005 ±0 


002 


0.007 


0.6 ±0.3 


2.23 


-0 


16 


0.03 


SI 


0.0053 


271 


0.023 ±0.001 


-0 


018 ±0.001 


-0.005 ±0 


002 


0.013 


0.3 ±0.4 


2.35 


-0 


20 


0.02 


S5 


0.0053 


548 


0.015 ±0.006 


-0 


011 ± 0.000 


-0.005 ±0 


004 


0.012 


0.1 ± 0.2 


2.39 


-0 


25 


0.02 


S6 


0.0054 


1063 


0.013 ±0.003 


-0 


007 ± 0.001 


-0.006 ± 


009 


0.010 


0.1 ± 0.2 


2.70 


-0 


32 


0.01 


11 


0.0112 


26 


0.064 ± 0.003 


-0 


060 ± 0.002 


-0.002 ±0 


001 


0.009 


1.1 ± 1.2 


2.01 


-0 


06 


0.06 


12 


0.0105 


55 


0.029 ±0.007 


-0 


027 ± 0.004 


-0.002 ± 


004 


0.007 


-0.0 ± 1.2 


9.11 


-0 


06 


0.03 



helicity flux across the equator was possible, playing thus a 
measurable role; see Table [T] 

To characterize the magnitude of the magnetic helicity, 
we give its value averaged over the range \z\ ^ L*. To com- 
pare this value with that from advective magnetic helicity 
fluxes, we should multiply the table entry for V ■ J-i by L» , 
which is about half the full vertical extent of the domain. 
Note that V • J-t and k\T z are actually comparable, even 



though T % can have no effect in the present geometry and 
gives zero divergence. 

We recall that j ■ b and a ■ b are approximately pro- 
portional to each other. This is also borne out by the 
present simulations where fc^ ff = j ■ b/a ■ b is constant and 
keff/k{ « 2. This confirms earlier findings of MCCTB and 
HB, where a similar value of k e ff was found. U nder isotropic 
cond itions, this ratio is approximately unity l|Brandenburd 
l200lf ). However, for Models N3 and N4, the correlation be- 
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Figure 5. Contributions to the magnetic helicity flux for 
Model W3. Upper panel: vertical profiles of magnetic helicity 
fluxes of the fluctuating field (solid line), compared with the con- 
tribution from the mean flow (dashed line). Lower panel: resid- 
ual between the two aforementioned fluxes (solid line) compared 
with a fit to the gradient of the magnetic helicity density from 
the small-scale field (dashed line). The fluxes are normalized by 
^"o = »7t0-B? q - 




0.0 L 



0.000 0.002 0.004 0.006 0.008 0.010 

s w 

Figure 6. Root-mcan-square value of the mean magnetic field, 
Brms, as function of S w for models N3, W2, M2, S2 and 12, which 
have R.m 50. 



tween j ■ b and a ■ b is poor, giving formally a negative value, 
so fceft is given as imaginary in Table [2] 

The quantity j ■ b/k[B^ q is systematically below unity, 
suggesting that the dynamo can only be expected to produce 
mean fields where B « B^ q . Finally, we also give the values 



of the flux divergence of the mean field V-^ m . These values 
are typically about 10 times larger than the flux divergence 
of magnetic helicity of the small-scale field, V • but it 
is of course only the latter that is relevant for alleviating 
catastrophic quenching. 

All simulations with wind show that the rms value of 
the mean field, B rms , declines slowly with increasing wind 
speed; see Figure [6] This result might just be a consequence 
of a gradual increase of the critical value of Rm above which 
dynamo action is possible. However, it could also be an in- 
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Figure 7. -B rms as a function of Rm in absence (solid line, 
Models Tl, T2, N1-N6) and in presence (dashed line, Models Sl- 
S6) of advective flux. The two dotted lines give the slopes —0.5, 
—0.4, and —0.17 for orientation. 
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Figure 8. Space-time diagrams of B y for wind Uo = 0.015e s 
(corresponding to Sw ~ 0.0055) for Models SI— S6 for different 
magnetic Reynolds numbers. From the top: Rm = 24, 51, 133, 
271, as well as 547 (bottom left) and 1063 (bottom right). 



dication that a fraction of the mean magnetic field is being 
removed from the d omain by the flow - as found in the 
mean- field models of IShukurov et all (|2006t ). 

In Figure [7] we see how B r ms decreases with increasing 
Rm. The scalings Rm" 0,4 and Rm" 0,17 are given for ori- 
entation and show that in the presence of advection -Brms 
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Figure 9. Scaling properties of the vertical slopes of 2£ ■ B, 
—2r]i±a j ■ b, and —V • for Models W1-W5 (upper panel) and 
for Models S1-S6 (lower panel). (Given that the three quanti- 
ties vary approximately linearly with 2, the three labels indicate 
their non-dimensional values at k\z = 1.) The second panel shows 
that for a stronger wind the contribution from the advective term 
becomes approximately independent of Rm for Rm > 170 (blue 
line), while that of the resistive term decreases approximately 
like Rm -2 ' 3 (red line), and 2£ ■ B decreases approximately like 
Rm- 1 / 2 (black line). 



varies much slower than Rm -1 , which is the slope antici- 
pated from catastrophic quenching m odels without a wind 
l|Brandenburg fc Subramanianl |2005|). Note, however, that 
DNS always gave a shallower slope (jBrandenburg fc Doblerl 
l200ll ) and, at larger values of Rm, B lms may have been a l- 
ready independent of Rm (|Hubbard fc Brandenburgll2012h . 
Indeed, without a wind (Sw = 0) the Rm dependence is 



compatible with a steeper Rm -1 ' 2 law, but it is less certain 
in this case. Looking at Figure [8l we can also see that there 
is no significant change of the cycle period with Rm. The 
high-resolution runs with Rm = 544 and 1061 are too short 
to cover a magnetic cycle, but one can see that the slope of 
the structure, which corresponds to the speed of the dynamo 
wave, is approximately unchanged. In the high-Rm models 
the fluctuations are more pronounced, but the peak-to-peak 
contrast is about the same for all runs. 

Table[5]shows that 2£B, 2r]fio j ■ b, and V-^f balance 
approximately to zero, confirming that the results represent 
a statistically steady state. All three quantities have approx- 
imately the same (nearly linear) 2 dependence for \z\ < £*, 
so that also the values of their three slopes must balance 
to zero, which is indeed the case. In Figure [9] we show the 
scaling properties of the aforementioned quantities for Mod- 
els W1-W5 and S1-S6. For Rm < Rm*, where Rm* w 170 
for Sw = 0.002 and Rm, « 120 for Sw = 0.005, the first 
two quantities decrease approximatively like Rm -1 , while 
the latter decreases only like Rm' 1 ' 2 , which is in agree- 



Table 3. Additional parameters of the simulations including Rm, 
magnetic diffusivity, the ratios 'B/u' (= B rms /it rms ) and '£?/&' 
(= Brms/trms), as well as Mach number and number of mesh 
points. 
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5.0 x 


10" 
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0.58 


1.13 





18 


64 


T2 


23 


2.0 x 


10" 


3 


0.48 


0.87 
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64 


Nl 


37 


1.0 x 


10" 


3 


0.53 


0.70 





15 


64 


N2 


81 


5.0 x 


10" 


-4 


0.58 


0.73 





16 
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N3 
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10" 


-4 


0.27 


0.33 





17 
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Nl 


397 


1.0 x 


10" 


-4 


0.27 


0.33 





16 


512 


N5 
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5.0 x 


10" 


-5 


0.18 


0.21 





15 


1024 


N6 


1073 


2.5 x 


10" 


-5 


0.11 


0.15 





11 


1024 


Wl 


24 


1.0 x 


10" 


-3 


0.61 


0.69 





10 


128 


W2 


51 


5.0 x 


10" 


-4 


0.42 


0.48 





10 


128 


W3 


129 


2.0 x 


10" 


-4 


0.36 


0.39 





10 


256 


W4 
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1.0 x 


10- 


-4 


0.28 


0.31 
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512 


W5 
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5.0 x 


10" 


-5 


0.19 


0.22 





11 


1024 
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51 


5.0 x 


10- 


-4 


0.36 


0.45 





10 


128 
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1.0 x 


10- 
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0.40 


0.55 





10 
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5.0 x 


10- 
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0.31 


0.42 
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128 
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10- 
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0.20 


0.27 
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256 
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1.0 x 


10- 
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0.17 


0.23 
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512 


S5 


548 
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10- 


-5 


0.15 


0.19 





11 


1024 


S6 
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10- 


-5 
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0.18 
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0.13 
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11 
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ment with the values ob tained by HB; see also Figure 10 of 
ICandelaresi et all (|201ll ) for a corresponding plot. 

However, for Rm > 200 the scaling of 2£ • B changes 
into an Rm -1 / 2 scaling; V ■ J-i is at first below 2r\\ioj ■ b, 
but for high enough Rm increases to reach an absolute value 
similar to that of 2£ B. This suggests that the simple expec- 
tation based on the naive extrapolation given from a linear 
fit is misleading, and that catastrophic quenching might be 
alleviated already for Rm > 1000. In the absence of a wind 
and for large magnetic Reynolds numbers (Models N4-N6), 
the divergence of the magnetic helicity flux shows strong 
fluctuations about zero, making it harder to determine an 
accurate magnetic helicity balance of small-scale fields. 

In Table [3] we summarize additional output parame- 
ters of the simulations including Rm, the magnetic diffusiv- 
ity, the ratios of the rms values of mean field to fluctuat- 
ing velocity and fluctuating magnetic field, i.e., B r ms/wrms 
and .Brms/fcrms, respectively, as well as Mach number and 
number of mesh points. As was already obvious from Fig- 
ure[7l -Brms/itrms (which is the /B cq ), decreases 

with increasing Rm, and the same is also true of the ratio 
Brms/^ms- The numerical resolution in the x direction, N x , 
is given in the last column. This is also the resolution used 
in the y direction, while that in the z direction is always 
twice as large. 



4 CONCLUSIONS 

In the present work we have examined the effects of an ad- 
vective magnetic helicity flux in DNS of a turbulent dy- 
namo. The present simulations without shear yield an os- 
cillatory large-scale field owing to the spatially varying ki- 
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netic helicity profile with respect to the equatorial plane. We 
emphasize in this context that t he possibility of oscillatory 
dynamos of a 2 type is not n ew (|Barvshnikova fe Shukurov! 
1 19871 : iRadler fc Brauerl fl987l ) . but until recently all known 
examples were restricted to spherical shell dynamos where a 
changes sign in the radial direction. The example found by 
iMitra et all l|2010bl ) applies to a spherical wedge with latitu- 
dinal variation of a changing sign about the equator. Similar 
results have also been obtained in a mean-field dynamo wit h 
a linear variation of a(z) cx z ([Brandenburg et al.ll2009T ). 
Our present simulations are probably the first DNS of such 
a dynamo in Cartesian geometry. Closest to our simulations 
are those of MCCTB who used perfectly conducting outer 
boundary conditions without wind, and also found oscilla- 
tory solutions. Surprisingly, however, oscillations are here 
only obtained if there is at least a slight outflow. 

One would have expected that catastrophic quenching 
can be alleviated if magnetic helicity is removed from the 
domain at a rate larger than its diffusion rate, that is, the 
advective term V • J-f dominates over the resistive term, 
2i)/xo j ■ b. Figure [9] shows that, for Rm < 200, the latter 
term decreases linearly with decreasing 77, while the for- 
mer only decreases proportional to r] 1 ^ 2 , i.e., proportional 
to Rm -1 / 2 . This would have led us to the estimate that for 
Rm w 4 • 10 3 the catastrophic quenching can be alleviated 
by a wind with Sw ~ 0.005. Our new results suggest that 
this can happen already for smaller values of Rm. The rea- 
son for this is still unclear. It is possible that catastrophic 
quenching was an artefact of intermedia t e valu es of Rm, 
as suggested bv lHubbard fc Brandenburd (|2012l '). or that a 
magnetic helicity flux can have an effect even though it is 
weak compared with diffusive terms. 

Finally, we should emphasize that we have only exam- 
ined here the case of subsonic advection. In real astrophys- 
ical cases, like galactic and stellar winds, the outflow is in- 
stead supersonic and can, thus, play an even more impor- 
tant role in alleviating the catastrophic quenching through 
the advection of magnetic helicity. This assumes, of course, 
that the dynamo is strong enough to be still excited in the 
presence of a stronger wind. 
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